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ABSTRACT 

We report on the first phase of our study of slightly rotating accretion 
flows onto black holes. We consider inviscid accretion flows with a spherically 
symmetric density distribution at the outer boundary, but with spherical 
symmetry broken by the introduction of a small, latitude-dependent angular 
momentum. We study accretion flows by means of numerical 2D, axisymmetric, 
hydrodynamical simulations. Our main result is that the properties of the 
accretion flow do not depend as much on the outer boundary conditions (i.e., the 
amount as well as distribution of the angular momentum) as on the geometry of 
the non-accreting matter. The material that has too much angular momentum 
to be accreted forms a thick torus near the equator. Consequently, the geometry 
of the polar region, where material is accreted (the funnel), and the mass 
accretion rate through it are constrained by the size and shape of the torus. Our 
results show one way in which the mass accretion rate of slightly rotating gas 
can be significantly reduced compared to the accretion of non-rotating gas (i.e., 
the Bondi rate), and set the stage for calculations that will take into account 
the transport of angular momentum and energy. 

Subject headings: accretion - hydrodynamics - black hole physics - outflows ~ 
galaxies: active - methods: numerical 
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1. Introduction 



Some of the most dramatic phenomena of astrophysics, such as quasars and powerful 
radio galaxies, are most likely powered by accretion onto supermassive black holes (SMBHs). 
Nevertheless, SMBHs appear to spend most of their time in a remarkably quiescent state. 
SMBHs are embedded in the relatively dense environments of galactic nuclei, and it is 
natural to suppose that the gravity due to an SMBH will draw in matter at high rates, 
leading to a high system luminosity. However, this simple prediction often fails as many 
systems are much dimmer than one would expect. 

To illustrate a key problem in constructing theoretical models for accretion onto a 
black hole, let us express the luminosity due to accretion as 

L ^ rjc^M,, (1) 

where we invoke the simplest assumption, that the luminosity is proportional to the mass 
accretion rate. Ma, and an efficiency factor, rj. The accretion luminosity is very uncertain 
because rj is uncertain: t] ranges from ~ 10~^ in a standard, radiatively efficient thin 
disk, to ~ 10~^^ for spherically symmetric accretion from a low density medium (e.g., 
Shakura & Sunyaev 1973; Shapiro 1973; Meszaros 1975). The mass accretion rate is also 
a source of uncertainty in estimates of the accretion luminosity because Ma depends on 
the physical conditions and geometry at large distances from the black hole. Nevertheless, 
it is customary to adopt the analytic formula due to Bondi (1952) to estimate the mass 
accretion rate. In his classic paper, Bondi (1952) considered spherically symmetric accretion 
from a non-rotating polytropic gas with uniform density p^o and sound speed Cqo at infinity. 
Under these assumptions, a steady state solution to the equations of mass and momentum 
conservation exists with a mass accretion rate of 

Mb = A47ri?|pooCoo, (2) 

where A is a dimensionlcss parameter that, for the Newtonian potential, depends only on 
the adiabatic index. The Bondi radius, Rb, is defined as 

Rs = (3) 

where G is the gravitational constant and M is the mass of the accretor. 

The Bondi accretion formula predicts that SMBHs in typical galaxies should be more 
luminous than observations indicate when rj is assumed to be as large as in a standard, 
radiatively efficient thin disk (e.g., Di Matteo et al. 1999, 2000, 2001; Loewenstein et al. 



-3- 



2001; Baganoff et al. 2001). In the context of equation (1), one possible explanation for 
this disagreement is that the black hole accretion flow can be radiatively inefficient because 
binding energy dissipated in the gas is advected through the event horizon before being 
radiated (Ichimaru 1977; Rees et al. 1982; Narayan & Yi 1994, 1995; Abramowicz et al. 
1995). However, pure advection-dominated inflow may not be the whole story. Even before 
recent observations forced us to confront very low SMBH luminosities, theorists had begun 
to realize that rotating, radiatively inefficient hydrodynamical (HD) flows are subject to 
strong convection (Begelman & Meier 1982; Narayan & Yi 1995), which can fundamentally 
change the flow pattern and its radiative properties (Blandford & Begelman 1999; Quataert 
& Narayan 1999; Narayan, Igumenshchev & Abramowicz 2000; Quataert & Gruzinov 2000). 
The theoretical studies showed that convection alters the steep (oc r~^/^) density proflles of 
the advection-dominated flow into a much flatter (oc r^^/^) proflle, which can explain the 
faintness of many SMBHs because it predicts relatively low density close to the black hole 
(i.e.. Ma is low in eq. 2). Similar structural changes occur in the magnctohydrodynamical 
(MHD) limit (Stone & Pringle 2001; Hawley, Balbus, & Stone 2001; Machida, Matsumoto 
& Mineshige 2001; Igumenshchev & Narayan 2002; Hawley & Balbus 2002), although 
here the turbulence is probably driven by magnetorotational instability (MRI) rather than 
thermal convection (Balbus & Hawley 2002; but see Abramowicz et al. 2002 and Narayan 
et al. 2002 for alternative views). 

The turbulent character of both HD and MHD models does not settle the issue of what 
happens to the energy and angular momentum that must be transported away. There are 
two possibilities: (i) turbulent transport effectively shuts off the accretion flow, turning it 
into a closed circulation (Narayan et al. 2000; Quataert & Gruzinov 2000) or (ii) turbulent 
transport drives powerful outflows that can strongly modify the black hole's environment 
(Narayan & Yi 1994, 1995; Blandford & Begelman 1999). Recent MHD simulations bring 
new insights that may help us to resolve this issue. For example, Hawley & Balbus's (2002) 
three-dimensional MHD simulations show that, with and without resistive heating, mass 
and energy in nonradiative accretion flows are carried off by an outflow in keeping with the 
outline of the second possibility. 

Another possible solution to the problem of very low SMBH luminosity is that mass 
is captured into the accretion flow at a rate that is far lower that Mb- Thus, a complete 
formulation of the accretion flow must also consider the outer boundary conditions. SMBHs 
draw matter from an extended medium and most authors assume that the Bondi (1952) 
formula provides an adequate approximation for the rate of mass supply. There has been 
little systematic work done to demonstrate that this assumption is justifled and correct. 
Igumenshchev & Narayan (2002) showed that even non-rotating Bondi accretion can be 
altered, in particular that the mass accretion rate can be reduced below the Bondi rate. The 
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cause of the Ma reduction in Igumenshchev & Narayan's simulations is the local generation 
of energy through by magnetic reconnection, which leads to the development of efficient 
convection. However, in Igumenshchev & Narayan's simulations the flow was supersonically 
injected into the computational grid through the outer boundary, at a specified rate 
(determined from the Bondi formula). Their calculations did not address self-consistently 
the problem of the mass accretion and supply rate because the latter was fixed at the outer 
boundary. As we mention above, the Bondi formula has been derived under the assumption 
that the gas is non-rotating and only under the influence of the central gravity. Thus, 
for a given gravitational field, the gas internal energy determines the accretion rate. By 
relaxing this assumption, introducing additional forces or sources of energy, one may find 
that the mass supply rate is much lower than the one predicted by the Bondi formula. For 
example, the rate at which matter is captured by a black hole can be severely limited when 
the matter is heated by X-rays produced near the black hole (Ostrikcr et al. 1976) or by 
mass outflow from the central region (Di Mattco ct al. 2002). In these two cases, the gas 
internal energy is increased. Introducing kinetic energy to the gas at infinity may have a 
similar effect: although the flow outside the Bondi accretion radius often can be described 
as nonrotating, even a tiny amount of angular momentum, I — when followed inward — 
could severely limit the rate at which matter is captured by the black hole. 

Our focus shall be on assessing the gross properties of rotating accretion flows onto 
black holes. We consider a classic Bondi accretion flow with the only modiflcations 
being the introduction of a small, latitude-dependent angular momentum at inflnity and 
a pseudo-Newtonian gravitational potential. The imposed angular momentum is small 
enough to have a negligible effect on the density distribution at the outer boundary, which 
remains spherically symmetric. Therefore, matter near the rotational axis can be accreted. 
We thus consider a very simple model of an accretion flow, far simpler than those occurring 
in nature, as we neglect the gravitational field due to the host galaxy, radiative heating and 
cooling, viscosity and MHD effects. For example, we will not consider here the transport 
of energy and angular momentum outward, as needed to accrete matter with a specific 
angular momentum higher than 2Rsc (where Rs — 2GM/c^ is the radius of a Schwarzschild 
black hole). Nevertheless, the results presented here provide a useful exploratory study of 
accretion onto black holes. In particular, our results constitute a "baseline" for evaluating 
the effects of dissipative and transport processes in subsequent work. They also serve as 
a "proof-of-conccpt" for the reduction of the mass accretion rate due to a small angular 
momentum in the accretion fiow. 



1.1. Expectations 
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Before we embark on a detailed analysis and numerical HD simulations, we first 
consider the problem of accretion of low-/ material in a general way. If the matter far 
from the SMBH (well outside the Bondi radius) has a uniform density and a specific 
angular momentum, /, which exceeds 2Rsc — a tiny value compared to the Keplerian 
angular momentum at Rb — then no accretion will take place without angular momentum 
transport. After a transient episode of infall, the gas will pile up outside the black hole 
and settle into a nearly steady state atmosphere bounded by a centrifugal barrier near 
the rotation axis. Reahstically, matter far from the SMBH will have a range of angular 
momenta, and in a steady state with axisymmetry, there will always be low-Z material close 
to the axis that can accrete steadily through a funnel along the rotational axis. 

For this highly idealized problem, one would expect that the mass accretion rate should 
scale with the dependence of the angular momentum on the polar angle, 6, at the outer 
radius. To'- the larger the amount of the material with I > 2Rsc at To, the lower the mass 
accretion rate. This geometrical argument on the nature of the Ma vs. I relation can be 
quantified as follows. If the angular momentum depends on 9 as 

m = iof{e), (4) 

where / = 1 on the equator (9 — 90°) and monotonically decreases to zero at the poles 
{9 = 0° and 180°), then a naive expectation would be that Ma/ Mb scales with the solid 
angle within which / < 2Rsc: 

^— = — — = l-cos6'o, (5) 
Mb ^t^ ^ ' 

where 9o is the width of the angular distribution for which I < 2Rsc. The latter can be 
formally defined as 



Oo = f 



-1 



mm 



(6) 



2Rsc\- 
k 

where /^^ represents the functional inverse of /. This simple geometrical argument is based 
on the assumption of radial flow and implies that if / < 2Rsc the material will be accreted 
approximately at the Bondi rate. If relation (6) were true then the accretion rate should 
decrease with increasing Iq for a fixed Rs- Additionally, one would expect that Ma/ Mb 
should be independent of Rb for fixed Rs and Iq. However, to determine the actual mass 
accretion rate even in this idealized case we need to perform numerical simulations, as 
just HD effects of the inviscid fiuid make the above geometrical argument invalid both 
quantitatively and qualitatively. 

The main result of our numerical HD calculations is that the properties of the accretion 
fiow do not depend as much on the outer boundary conditions (i.e., the amount as well as 
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distribution of the angular momentum) as on the geometry of the non-accreting matter. 
Material with / 2Rsc cannot accrete and forms a thick torus near the equator. This 
thick torus and its formation have been a subject of numerous studies (see below). Our 
simulations show that the dependence of angular momentum on 6 in the torus gets weaker 
with decreasing radius. The material with I ^ 2i?5C inflows in the polar region, turns 
around as it reaches a centrifugal barrier, and then starts to outflow along the equator. 
As a result, a thick torus of nearly uniform specific angular momentum forms as gas of 
the highest angular momentum {I Iq > 2Rsc) near the equator is replaced by gas of 
lower angular momentum {I ^ 2Rsc). The geometry of the polar region, where material is 
accreted (the funnel) and the mass accretion rate through it are constrained by the size and 
shape of the torus. We describe the size of the torus by an angle, 9t, between the torus's 
upper envelope and the pole at a characteristic radius. 

Our HD models show that the vs. I relation has three regimes for a given f{d): (i) 
for low I (i.e., I < 2Rsc for all 9 at large radii), the torus does not form and Ma = Mb, 
(ii) for intermediate /, or more appropriately where there is a narrow range of 9 for which 
/ > 2RsC so 9o > 9t, Ma ~ const with the actual value of the constant depending on the 
ratio Rs/Rb and (iii) for high /, or in the case for which / > 2Rsc at nearly all 9 so 9o < 9t, 
Ma decreases with increasing /q. 

Comparing the Ma vs. I relation based on our HD models with that described by 
equation (5), we find that the two relations agree exactly in the first regime, disagree 
qualitatively and quantitatively in the second regime [the HD models predict Ma lower 
than eq. (5)], and agree again but only qualitatively in the third region [the HD models 
predict Ma higher than eq. (5)]. Thus the geometrical argument used above does not 
hold. However, we can use a modified version to describe the key aspects of the vs. I 
relation. The modification to the geometrical argument is to replace the solid angle within 
which I < 2Rsc at the outer boundary, AQq, with the solid angle within which I < 2Rsc 
at a characteristic radius comparable with the sonic radius, AQf (i.e., the solid angle of 
the accretion funnel). Thus, Ma/ Mb ~ AQf /An. The insensitivity of Ma to the angular 
momentum at infinity, in the second regime, can be attributed to the relative insensitivity 
of the torus, and thus the funnel, to the angular momentum distribution at infinity. In 
terms of the sohd angle within which / < 2Rsc, this corresponds to Aflf = constant for 
variable Iq, provided that AQf < AQg (i.e., 9o > 9t). On the other hand, the decrease of 
Ma with increasing Iq, in the third regime, can be attributed to the fact that AQf decreases 
with increasing Iq, provided that AQf > AQ^, {9o < 9t). We find that the mass accretion 
rate decreases with increasing Iq more slowly than predicted by eq. (5) in the third regime 
because the sonic point radius starts to increase as the funnel gets narrower. 
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1.2. Previous Work 

Similar calculations have been performed before. For example, the formation of 
rotationally supported thick tori from inviscid accretion of gas with various initial angular 
momentum distributions has been reported (Hawley, Smarr & Wilson 1984a; 1984b; Clarke, 
Karpik & Henriksen 1985; Hawley 1986; Molteni, Lanzafame & Chakrabarti 1994; Ryu et 
al. 1995; Chen et al. 1997). However, there is one key difference between our work and 
some past work: our outer radial boundary is located outside the Bondi radius and we adopt 
subsonic, Bondi-like outer radial conditions whereas Molteni et al. 1994, Ryu et al. 1995, 
and Chen et al. 1997 (see also Toropin et al. 1999; Kryukov et al. 2000; and Igumenshchev 
& Narayan 2002) imposed outer boundary conditions inside the Bondi radius or even inside 
the sonic radius. The latter approach allows one to study cases where Rs/Rb is as low as 
in some real systems (e.g., Rs/Rb — 10~^ in Chen et al. 1997) but this approach is not 
suitable for addressing our main issue: what is the mass supply rate. The approach adopted 
by Hawley et al. (1984a, 1984b), Clarke et al. (1985) and Hawley (1986) is much closer 
to ours as they also used subsonic outer boundary conditions. However, these authors did 
not consider how the accretion rate onto the black hole depends on the angular momentum 
distribution beyond the Bondi radius but rather focused on the formation of the thick 
torus. As far as the treatment of the outer radial boundary is concerned, our simulations 
are also similar for those of Ruffert (1994), who studied three-dimensional hydrodynamic 
Bondi-Hoyle accretion. 

Other studies are also relevant to our work. Several authors considered accretion onto 
black holes with a focus on the evolution of rotationally supported thick tori including 
the transport of angular momentum and energy (e.g., Igumenshchev & Abramowicz 1999; 
Stone, Pringle & Begelman 1999; Stone & Pringle 2001; Machida et al. 2001; Hawley & 
Balbus 2002). The main difference between our work and these studies is that the other 
authors adopt the point of view that virtually all of the material at large radii has too much 
angular momentum to be accreted without the transport of angular momentum. They 
assume that material with zero or very low angular momentum is unimportant dynamically, 
and that accretion is dominated by angular momentum and energy transport processes. 
For example, for their initial conditions Stone et al. (1999) and Hawley & Balbus (2002) 
adopted a bounded torus in hydrostatic equilibrium with constant angular momentum, 
embedded in zero angular momentum ambient gas which is also in hydrostatic equilibrium. 
Thus, these calculations were set up so that, if not for the transport of angular momentum, 
there would be neither time evolution nor mass accretion. We recognize that transport 
processes are essential, but assert that the zero or very low angular momentum material 
also deserves a rigorous treatment because it can play an important role in determining the 
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total mass supply and accretion rate. 

In this paper, we consider a far simpler case of an accretion flow than those occurring 
in nature (see Section 4). For example, we neglect viscosity and MHD effects. In particular, 
the MRI has been shown to be a very robust and universal mechanism to produce turbulence 
and the transport of angular momentum in disks at all radii (Balbus & Hawley 1998). 

The outline of this paper is as follows. We describe our calculations in Section 2. In 
Section 3, we present our results. We summarize our results and discuss them together with 
their limitations in Section 4. 



2. Method 



2.1. Hydrodynamics 




To calculate the structure and evolution of an accreting flow, we solve the equations of 
hydrodynamics 

^ + PV.V = 0, (7) 
p^ = -VP + pV$, (8) 

= -PV • V, (9) 

where p is the mass density, P is the gas pressure, v is the velocity, and e is the internal 
energy density. We adopt an adiabatic equation of state P — (7 — l)e, and consider models 
with 7 = 5/3. Our calculations are performed in spherical polar coordinates {r,9,(f)). We 
assume axial symmetry about the rotational axis of the accretion flow {9 — 0°). 

We present simulations using the pseudo-Newtonian potential $ introduced by 
Paczyhski & Wiita (1980) 

* = -^. (10) 
r-Rs ^ ' 

This potential approximates general relativistic effects in the inner regions, for a nonrotating 

black hole. In particular, the Paczynski-Wiita potential reproduces the last stable circular 

orbit at r = 3Rs as well as the marginally bound orbit at r = 2Rs. 
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2.2. Initial conditions and boundary conditions 



For the initial conditions we adopt a Bondi accretion flow with zero angular momentum 
everywhere except for the outermost part of the flow. In what follows we briefly review 
the basics of Bondi accretion that allow us to specify details of our initial and boundary 
conditions as well as to interpret our results. 



The Bernoulli function can be written as 



+ 



/2 



2r2 sin^ 9 



(11) 



where H is the enthalpy and I = v^r sin 6 is the specific angular momentum. For a polytropic 
equation of state P — Kp'^, the polytropic constant can been expressed as K — p^^^c^/j, 
where Coo is the sound speed at infinity [i.e., = {dP/dp)oo\- Therefore the enthalpy 
becomes 



H 



-I 



p dP 

Poo P 



7 



7 



7 



.Poo, 



7-1 



(12) 



Let us consider spherically symmetric, steady-state Bondi accretion onto an object 
with a Paczyhski-Wiita (PW, hereafter) potential. In such a case, ve — I — ^ and the 
Bernoulli function simplifies to 



B = 



GM 



+ 



Rs 7 



.Poo, 



>-l 



- 1 



+ 



0. 



(13) 



Introducing dimensionless variables a = p/poo, x = r/Rs, R's = Rs/Rb, and 
Ur = Vr/coo, wc rcwritc eq. (13) as 

1 1 ?/2 

B' = —- + ^(a^-i _ 1) + = 0. 

x-R's 7-r ' 2 



(14) 



The transonic solution of the Bondi problem yields the sonic point at 

5 - 37 + 8R's + 



Xc 



(15) 
(16) 



where 

A = (5 - 37 + 8R'sf - 64(1 - 7 + R's)R's- 

The dimensionless mass accretion rate A = x^aur for the transonic solution can be 
calculated as 

X ^ xla,Ur,s ^ xla(;+'y\ (17) 
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with the dimensionless density at the sonic point, Ug, given by 




(18) 



For 7 = 5/3, the dimensionless sonic radius and mass accretion rate are 



xg^R's + ^/wys 



(19) 



and 




(20) 



respectively. Note that, for the PW potential as well as the fully general relativistic problem 
(Bcgclman 1978), Bondi flow with 7 = 5/3 has a sonic point at roughly the geometrical 
mean between the Bondi radius and the Schwarzschild radius. This is in contrast to the 
purely Newtonian case in which Xs = {5 — 37) /4 — >• as 7 — >• 5/3. 

We will use the results for Bondi accretion with the PW potential [i.e., equations 
(17)-(20)] as a reference point for analyzing our results and to set the initial conditions. 
Namely, we adopt vg = while Vj. and p are computed using the Bernoulli function and mass 
accretion rate Mb = —AiTr'^pVj. = XAnR'^p^Coo- We set p^o = 1 and specify Cqo through R'^ 
(note that R'g = 2c^/c^). We complete specification of the initial conditions by adopting a 
non-zero specific angular momentum Z for r > XsRb- We also ran a few models in which the 
specific angular momentum initially is non-zero only at the outer radial boundary Tq. We 
found that, except for the initial transient, models with these initial conditions for I give 
the same results. However, adopting a non-zero specific angular momentum / for r > XsRb 
reduces the computational time significantly because the initial transient lasts for a much 
shorter time. 

We consider a general case where the angular momentum depends on the polar angle 

via 



with / = 1 on the equator {9 = 90°) and / = at the poles {9 = 0° and 180°). We express 
the angular momentum on the equator as 



where R'^ is the "circularization radius" on the equator in units oi Rb for the Newtonian 
potential (i.e., GM/r"^ = vl/r at r = R'cRb)- 

We adopt three forms for the function f{9): 



l{ro,9)^lof{9), 



(21) 




(22) 



fi{9) = 1 - |cose 



(23) 
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hie) = 1 - cos 



,10 



(24) 



and 



h{0) 



for e < Oo and 6 > 180° - 
lo for 9o< 9 < 180° - 9o. 



(25) 



Our choice of the 6 dependence of I at infinity is motivated by the following: the I 
distribution described by eq. (23) allows us to study a case in which the flow at large radii 
is very similar to the Bondi flow as the material with / > 2RsC is conflned to within a 
narrow range of 9 above and below the equator if Iq is close to 2RsC. In such a case most 
of the material at large radii has / < 2Rsc and can be accreted, in principle. On the other 
hand, the I distribution described by eqs. (24) or (25) allows us to study the opposite 
case, in which the material with I > 2Rsc at large radii occupies nearly the entire range 
of 9 above and below the equator, even if Iq is not much larger than 2Rsc. Then most of 
the material at large radii has / > 2RsC] only in a narrow polar region is / low enough for 
accretion to take place {I < 2Rsc). We note that for the 9 dependence of / described by fi 
and /2, we calculate 9q using eq. (6) and an assumed value for Iq. For the 9 dependence of I 
described by /a, Iq as well as 6*0 are free parameters. 

Our standard computational domain is defined to occupy the radial range 
r,i = 1.5 Rs < T < To = 1.2 Rb and the angular range 0° < ^ < 180°. We consider 
models with R'^ from 10~^ to 10^^'^. The r — 9 domain is discretized into zones. For 
R'g = 10~^, 10~^, and 10~^'^, our numerical resolution consists of 100, 140, and 160 zones 
in the r direction, respectively. In the 9 direction, our numerical resolution consists of 100 
zones for all values of R'g. We usually fix zone size ratios, drk+i/drk = 1.05, d9i/d9i+i = 1.05 
for 0° < ^ < 90° (i.e., the zone spacing is decreasing in this region) and d9i+i/d9i — 1.05 for 
90° < ^ < 180°. For runs with very high Iq we adopt d9i/d9i+i = 1.0 for 0° < ^ < 180° 
(see section 3.3). For runs with a very narrow funnel in which we use a step function for 
f{9), we adopt d9i+i/d9i = 1.03 for 0° < ^ < 90° (i.e., the zone spacing is increasing in this 
region) and d9i/d9i+i = 1.03 for 90° < ^ < 180°. 

The boundary conditions are specified as follows. At the poles, (i.e., ^ = 0° and 
180°), we apply an axis-of-symmetry boundary condition. For the inner and outer radial 
boundaries, we apply an outflow boundary condition. Our choice for the location of the 
inner radial boundary, = 1.5Rs, ensures that the flow near this boundary is supersonic 
and the outflow condition is appropriate. To represent steady conditions at the outer radial 
boundary, during the evolution of each model we continue to apply the constraints that in 
the last zone in the radial direction, ve = 0, = lof{9)/rsm9, and the density is fixed 
at the Bondi value at all times. Note that we allow Vr to fioat. We have found that this 
technique, when applied to calculations of spherically symmetric accretion, produces a 
solution that relaxes to the steady-state Bondi solution for I — 0. For non-zero I, we find 
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that the outer radial boundary must be sufficiently far from the circularization radius to 
ensure that the assumption of pure radial flow (i.e., V0{ro) — 0) does not change the solution 
(see Section 3). 

To solve eqs. (7)- (9) we use the ZEUS-2D code described by Stone & Norman (1992), 
modified to implement the PW potential. As with all Eulerian codes, one ought to test how 
much numerical diffusion error ("numerical viscosity") affects the results. In particular, 
shocks are captured with a standard quadratic artificial viscosity (see Stone & Norman 
1992). Since we compute the evolution of adiabatic flows without any physical viscosity, 
numerical effects can, in principal, limit the accuracy of the results. The artificial viscosity, 
which is not shown in our equations (7)-(9), is present only when the (inward) velocity 
divergence is nonvanishing. When present, artificial viscosity heats the gas and transports 
angular momentum. As a test, we reran model A04fila with a resolution of 150 in the r 
direction and of 200 in the 6 direction (i.e., the resolution of the test run is twice that in 
model A04fila in both directions). We found that the mass accretion rate in model A04fila 
and its higher resolution counterpart agree to within one part in 10^. As for heating of 
the gas due to artificial viscosity, we found that the polytropic constant, K, is conserved 
in all of our simulations, except for departure of ^ 1% near the equator for small radii, 
where weak shocks appear. We conclude that our numerical simulations treat inviscid fiows 
satisfactorily. 



3. Results 

Our numerical models are specified by several parameters. The length scale is 
determined by the black hole radius in units of the Bondi radius, R'g. Our second parameter 
is the adiabatic index, 7. The last parameter (or a function rather) is the angular 
momentum at the outer radial boundary, I — lof{d)- 

In real systems, R'g is relatively small (10~^ ^ R's — 2c^/c^ ^ 10~^) because the 
sound speed at large distances from a black hole is very small compared to the speed of 
light. Models with so great a radial domain would be very demanding computationally. We 

therefore decided to explore the main features of our model by considering a smaller radial 
domain, lO^'^'^ < R'g < 10~^. We consider adiabatic fiows with 7 = 5/3. Finally, we assume 
an angular momentum distribution at the outer radial boundary as described in Section 2. 
We focus our attention on accretion of matter with low angular momentum, i.e., where the 
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corresponding centrifugal force is small compared to gravity for all 6 at the Bondi radius: 

'^/^W <,^^. (26) 



In the limit Rb ^ Rs and using our definition of R'q, we can rewrite the above equation as 

^ < -L (27) 

Table 1 summarizes the properties of the simulations we discuss here. Columns (2) 
through (7) give the numerical resolution in the radial direction; the black hole radius 
compared to the Bondi radius, R'g] the circularization radius compared to the Bondi radius, 
R'fj] the specific angular momentum on the equator at r = Tq, Iq, in units of 2Rsc; the 
width of the angular distribution for which I < 2Rsc, Bo] and the angular momentum 
dependence on the polar angle at the outer boundary, f{9), respectively. Table 1 also 
presents the final time at which we stopped each simulation (all times here are in units of 
the Keplerian orbital time at r = Rb), the maximum specific angular momentum at the 
inner radial boundary, Z™""^, and the time- averaged mass accretion rate through the inner 
radial boundary measured near the end of the simulation, in units of the corresponding 
Bondi accretion rate. Finally, column (11) gives comments about runs different from the 
standard runs (e.g., the outer radial boundary set at 12Rb instead of 1.2Rb) 

Our simulations show that for Zq > 2RsC the accretion flow consists of a thick, 
equatorial torus where the material has too much angular momentum to be accreted, and a 
supersonic polar funnel where the material has I low enough to be accreted. We describe 
an example of such an accretion flow in some detail flrst (Section 3.1). This is followed by a 
limited parameter survey in which we focus on varying three key aspects of our models: the 
maximum angular momentum for a given distribution, the angular distribution of angular 
momentum on the outer boundary, and the black hole radius compared to the Bondi radius. 



3.1. Accretion flow consisting of a torus and a funnel 



In this section we describe the properties and behavior of our model in which 
R's = 10"^, Rq = 10"^ and f{9) = 1 - | cos^| (model B04fla). For the above parameters, 
the speciflc angular momentum on the equator at the outer boundary is ~ 3.5 in units of 
2Rsc, thus the material that cannot be accreted onto the black hole is located relatively 
close to the equator, 45° < 9 < 135° at the outer boundary. 
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Figure 1 presents a sequence of specific angular momentum contours and velocity 
fields for model B04fla. The length of the arrows in the lower panels is proportional to 
the poloidal velocity \Jv^ + Vq. As we mentioned in Section 1, after a transient episode of 
infall, the gas with / > 2Rsc piles up outside the black hole and settles into a thick torus 
bounded by a centrifugal barrier near the rotation axis. The low-/ material close to the axis 
accretes almost steadily through a funnel. The distribution of specific angular momentum 
in the torus becomes increasingly homogeneous as material with I ~ 0.9 replaces material 
with higher I. The fiow in the torus is subsonic, variable and is directed outward near the 
equator and inward close to the poles. There is no symmetry with respect to the equator. 
The time-dependence persists even after Thus, a meridional circulation is established 

in the torus. However, it is important to note that the time-averaged gross properties of the 
flow (such as the mass accretion rate and the shape of the torus and funnel) settle down to 
a steady state. 

To show the accretion fiow in more detail. Figure 2 presents the enlargement of 
two panels from Figure 1 (the panels from the the end of our simulations). This figure 
shows also the radial sonic surface, that is, the location where the radial Mach number 
{M.f = Vr/cs) equals one. Note that the specific angular momentum in the torus for small 
radii is nearly constant (/ ~ 1 in units of 2Rsc). Additionally, for small radii the sonic 
surface approximately coincides with the / a; 0.9 x 2Rsc surface. The 'south' lobe of the 
sonic surface is more elongated than the 'north' lobe at the end of the simulations. However, 
in a time-averaged sense both lobes are very similar. 

Next we consider the angular dependence of the fiow at small radii. Figure 3 is a plot 
of the angular dependence of the Mach number, specific angular momentum and mass fiux 
density at r = = l.BRs- This figure shows three important features of the accretion fiow: 
(i) the accretion fiow is highly supersonic in the radial direction, (ii) Z^"^ is not 2Rsc but 
rather somewhat smaller (0.9 x 2Rsc) and (iii) the mass fiux density is a strong function of 
the polar angle: it is nearly fiat near the poles and peaks near the equator where I ~ 0.6 in 
units of 2Rsc. 

To provide more insight into the character of the accretion fiow. Figure 4 shows the 
angular dependence of fiow properties on the sonic surface (the Mr = 1 surface). The top 
panel shows the radius at which Mr = 1, r^. The second panel from the top shows the 
specific angular momentum and the latitudinal Mach number, Mg = Vg/Cg. The second 
panel from the bottom and the bottom panel show, respectively, the mass fiux density (pVr) 
and a measure of how much the fiow deviates from a purely radial fiow at the sonic surface, 
A^. We define AO as the difference between a given polar angle, 9, on the sonic surface, and 
the polar angle at which the gas with a given I originated at r = Tq. Because our models 
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conserve specific angular momentum along streamlines, can be formally estimated as 

M{d) = d-r\l{rs,e)/h), (28) 

for 0° < ^ < 90° and 

A^(^) = 180° - e - r\l{rs, e)/lo), (29) 

for 90° < ^ < 180°. 

Figure 4 shows that the sonic radius is the highest along the poles, comparable to the 
Bondi sonic radius [see eq. (19)] and decreases by approximately an order of magnitude 
near the equator (see also Figure 2). Additionally, Figure 4 shows that the flow deviates 
from radial flow on the sonic surface. The deviation from radial flow increases with 
increasing angle from the poles. For angles from the poles ^ 30°, the deviation is small and 
the radial approximation holds. However, the mass flux density is dominated by matter 
originating at angles from > 40° from the poles, for which the deviation from radial flow 
can be large (i.e., 10° < A^^ < 50°). We note that the total mass accretion rate is dominated 
by streamlines which deviate only moderately from radial flow on the sonic surface. This 
conclusion is based on the observation that A^ ^ 10° for material with / ~ 0.6 x 2Rsc. 

Finally, to provide some insight into the time dependence. Figure 5 shows the time 
evolution of the mass accretion rate in units of the corresponding Bondi rate. Initially, Ma 
drops from 1 to 0.23 at t = 0.25. Then it starts oscillating around 0.3 with an amplitude of 

~ 0.1 Ma/Ms. 



3.2. Dependence of accretion flow properties on Iq 

To check the trends of the accretion flow with the main model parameters, we have 
performed a limited parameter survey. We first describe our results for models with all 
parameters the same as in model B04fla except for Iq. 

The main result from our runs with various Iq is that the mass accretion rate does not 
depend on this parameter (see Table 1), within the constraints imposed by assuming low 
angular momentum at Rb [eq. (27)]. The basic reason for this result is the fact that the 
properties of the accretion flow are determined by the geometry of the torus, which has a 
uniform specific angular momentum slightly larger than the maximum I of the accreted 
gas. The material with / ^ 2i?5C inflows in the polar region, turns around as it reaches 
a centrifugal barrier, and then starts to outflow along the equator. As a result, a thick 
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torus forms with gas of the highest angular momentum (Z Zq > 2i?5c) near the equator 
being replaced by gas of lower angular momentum {I ^ 2Rsc). The geometry of the polar 
region, where material is accreted (the funnel) and the mass accretion rate through it are 
constrained by the shape of the torus. 

To illustrate the insensitivity of the inner flow to the angular momentum at large radii, 
Figure 6 shows the angular momentum contours for our model BOlfla at tf = 3.5. For 
this model we assumed Iq = 2RsC so that in the purely radial case all material should be 
accreted. However, the flow converges toward the equator, and gas pressure effects allow 
only the material with Iq ^ 0.87 x 2RsC to be accreted. The material with higher / turns 
around and forms a thick torus with properties very similar to those of model B04fla, which 
has a much higher I at large radii. Comparing the left-hand panel of Figure 6 with the top 
panels of Figure 1 shows clearly how different the angular momentum distributions are at 
large radii in these two runs. However, comparing the right-hand panel of Figure 6 with the 
left-hand panel of Figure 2 shows that the angular momentum distribution at small radii, 
the shape of the torus (the I — 0.9 contour) as well as the sonic surface are qualitatively 
similar for both runs. 

We conclude that the relative insensitivity of the torus to the angular momentum 
distribution at infinity is responsible for the fact that Ma does not depend on Zq and that 
the geometrical argument we used in Section 1 is invalid. As runs B04fla and BOlfla 
illustrate, the width of the / distribution for / < 2RsC at infinity does not control the mass 
accretion rate. Even if material with low / comes from close to the equator it still has to go 
'around' the inner torus. We suggest that a more useful way of geometrically determining 
the mass accretion rate is to compare the angle between the torus upper envelope and the 
pole, 9t, and the width of the angular distribution for which I ~ 2Rsc, 9o- However, such 
a method of determining Ma poses the difficulty that 9t depends on radius and we need to 
know at what radius we have to measure Of. 



3.3. Dependence of accretion flow properties on the angular distribution of I 

Motivated by the conclusion from the previous section, we have performed a few 
simulations using an angular distribution for I that yields a smaller value of 9o than that 
given by /i (sec cq. 23), for a given Iq. We wish to check whether the insensitivity of Ma to 
Iq, as described above, holds in cases where 9o is small, in particular when 9o < 9t. 

Inspection of the expression (27) shows that there are three possible ways of decreasing 
9o within the constraints imposed by the assumption of low angular momentum: (i) reduce 
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R'q by increasing Rb for given Rc and Rg and increase Iq at the same time, (ii) choose a 
function f{9) that increases with increasing 9 very quickly for smaU 9 and very sfowly for 
larger 9, or (iii) both. We will be able to explore the first possibility only in a limited way 
(see next section) because, as we mentioned above, models with Rb very high compared to 
Rs are very demanding computationally. However, we can explore the second possibility at 
relatively low computational cost for moderate Rb- 

As before, we consider in detail only one example: for our model B08f2a, all parameters 
are the same as for model B04fla with the exceptions that / scales with at To as 1 — cos^° 9 
and d9i/d9i^i = 1. Such a flat distribution of / for high 9 reduces the size of the polar 
region where I < 2Rsc (i.e., the 9o angle). Thus, we expect that 9o < 9t in this case and 
the mass accretion rate will be lower than in model B04fla. Indeed, we find that Ma in 
model B08f2a is lower than that for model B04fla by a factor of ~ 3. We observe a similar 
reduction of Ma for all runs using /2 compared to runs using /i, regardless of Rs/Rb- 

Our experiments with various Iq and angular distributions of / show that 9o appears 
to be the key parameter that determines the effect of the angular momentum at the outer 
boundary on the mass accretion rate through the inner boundary. To illustrate the role 
of 9o, Figure 7 plots the mass accretion rate as a function of 9o for all our runs with 
Rs/Rb = 10~^, without making any distinction between models with different f{9) and/or 
Iq. Figure 7 shows one of the key results of this paper. For large 9o — corresponding to a 
broad funnel at Rb — Ma is constant, whereas for small 9o — corresponding to a narrow 
funnel at Rb — Ma decreases with decreasing 9o- 

Figure 7 also compares our numerical results (the solid line) with the simple prediction, 
eq. (5) (the dashed line). Equation (5) predicts that Ma decreases monotonically with 
decreasing 9o for all 9. Generally, our prediction from Section 1 overestimates Ma for large 
9o and underestimates Ma, for small 9o. The decrease of Ma with decreasing 9o shown in our 
numerical simulations is consistent with the theoretical prediction, but only qualitatively 
and only in the limit of a narrow funnel. 

We note that in the limit of a narrow funnel the sonic radius increases with decreasing 
9o- We can understand the Ma vs. 9o relation in the limit of a narrow funnel by combining 
this increase in the sonic radius with two other facts: (i) the sofid angle subtended by the 
accretion funnel decreases with increasing radius and (ii) the solid angle subtended by 
the accretion funnel is greater than AQp (AQj = AQq at by definition). Thus, AQj 
and consequently Ma decrease with decreasing 9o. This explains why the decrease of Ma 
with decreasing 9o shown in our numerical simulations is qualitatively consistent with the 
theoretical prediction. The explanation for the quantitative difference is simply related to 
the second fact listed above, i.e., AOj > Af^o at the sonic radius, < Tq. 



-18- 



3.4. Dependence of accretion flow properties on Rs/Rb and To 

The most difficult parameters of our models to explore are the Rs/Rb ratio and the 
outer radius of the computational domain. As mentioned earlier, we have to consider values 
of Rs/Rb orders magnitude larger than those found in astrophysical objects. Similarly, 
exploring a wide range of is computationally prohibitive. Therefore we have performed 
only a few simulations in which we vary these parameters. We check then whether our 
results allow us to extrapolate an M vs. Rs/Rb relation to very low values of Rs/Rb- 

Our results from the runs with various To allow us to make a few observations. First, a 
proper study of Bondi accretion with low angular momentum requires Rs/Rb < 10^^. For 
Rs/Rb = 10^^ or larger, the assumption of low angular momentum (eq. 27) is valid only 
for a very narrow range of /q while the assumption of radial flow at the outer boundary is 
valid only for To much larger than the Bondi radius. For example, our models A03fla and 
A03flb show that the location of the outer boundary can change the solution dramatically 
for Rs/Rb < 10~^ and moderate Iq of 2.5 x 2Rsc. For these parameters the circularization 
radius is relatively high, R'q — 0.5, and the assumption that the flow is radial close to the 
Bondi radius is not valid. Our models A03fla and AOSflb differ only in the location of 
the outer radial boundary, = 1.2Rb and 12 Rb, respectively. However, this difference 
produces dramatic change in the accretion flow. For model AOSflb, the flow is non-radial 
even as far as the Bondi radius and material with / as low as ~ 0.4 x 2Rsc turns around 
near the equator and never makes it through the inner boundary. Since the point where this 
material turns is relatively far from the inner boundary, the sonic surface of the flow coming 
from the polar region is approximately spherical and almost as large as the sonic surface 
in Bondi accretion with I — 0. Because of these properties of the sonic surface the mass 
accretion rate is almost as high as the Bondi rate\ This model once again demonstrates 
that the mass accretion rate is determined more by the shape of the flow that does not 
accrete than by the angular distribution of / at large radii. On the other hand, model 
AOSfla has = 1.2 and forces the flow to be radial near the Bondi radius via our outer 
boundary conditions. Therefore, the flow converges toward the equator at smaller radii and 
the material that goes through the inner boundary has I higher than in the model A03flb. 
Additionally, the point where the material with too high I turns around near the equator is 
closer to the black hole than the sonic point in the polar region. Consequently, for model 
AOSfla the sonic surface has the 'figure-8' shape and the mass accretion rate is lower than 
the Bondi rate. 

Because of the sensitivity of the results to the location of the outer boundary for 
models with Rs/Rb — 10~^, we do not consider these models in great detail. However, we 
note that these models (with fixed Tq) show an vs. I relation very similar to the one we 
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found for models with Rs/Rb — ^0 ^. 

We base most of our observations and conclusions on models with Rs/Rb ~ 10~^ not 
only because they are computationally affordable but also because they do not depend on 
the location of the outer boundary. Comparing models B04fla and B04flb we find that 
moving To from 1.2 Rb to 12 Rb does not change the solution for the accretion flow. The 
reason for this is simply that, contrary to models with Rs/Rb = 10~^, the Bondi radius 
is large enough compared to Rc that the 6 component of centrifugal force is relatively 
small. Consequently, the radial approximation at the Bondi radius is valid for models with 
Rs/Rb = 10"^ 

Unfortunately, models with the ratio Rs/Rb lower than 10~^ are computationally 
too expensive to be explored as fully as models with Rs/Rb = 10~^. For example, it 
takes about three weeks of CPU time on a modern workstation to run one model with 
Rs/Rb = 10~^'^ for 3 dynamical time scales at the Bondi radius. Therefore, we limited 
ourselves to two models with 10~^'^, model C04fla and C08f2a. Our primary goal of 
running these models was to check whether the mass accretion rate would change compared 
to models with Rs/Rb = 10~^. 

Model C04fla is meant to represent a model with a wide accretion funnel (the second 
regime in the Ma vs. I relation) while model C08f2a is meant to represent a model with 
a narrow accretion funnel (the third regime in the vs. I relation). For model C04fla, 
the mass accretion rate is smaller than the corresponding model with Rs/Rb — 10~^ (i.e., 
model B04fla, see Table 1) but we note that the mass accretion rate was increasing with 
time at the end of our calculations. Judging from the trend of Ma with time we conclude 
that the mass accretion rate docs not change much when the ratio Rs/Rb decreases from 
10^'^ to lO^'^'^ in the wide funnel case. We reach the same conclusion for the models in 
the narrow funnel case. Wc recognize that our conclusions arc based on very limited data 
but they are consistent, at least for small Rs/Rb, with our understanding of how the mass 
accretion rate is determined. 



4. Discussion 

This paper presents the first phase of our study of shghtly rotating accretion flows 
onto black holes. We decided to consider a far simpler case of an accretion flow than 
those occurring in nature. In particular, we neglected the gravitational field due to 
the host galaxy, radiative heating and cooling, viscosity and MHD effects. Perhaps the 
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most important simplification we made is neglecting the transport of energy and angular 
momentum outward as needed to accrete matter with a specific angular momentum higher 
than 2Rsc. Nevertheless, our results provide a useful exploratory study of accretion onto 
black holes as they have revealed unexpected properties and complexity of accretion flows 
in even this simple case. Clearly, a lot more work is needed to give a definitive answer to 
the question of whether slow rotation of gas at large radii is really enough to reduce the 
mass accretion rate to the level required by observations. In what follows we will summarize 
our results, briefly review the limitations of our work, and discuss how the physical effects 
neglected here may change the results. 

We have performed numerical 2D, axisymmetric, hydro dynamical simulations of 
slightly rotating, inviscid accretion flows onto a black hole. We attempt to mimic the 
boundary conditions of classic Bondi accretion flows with the only modifications being the 
introduction of a small, latitude-dependent angular momentum at the outer boundary and 
a pseudo-Newtonian gravitational potential. The distribution of I with latitude allows the 
density distribution at infinity to approach spherical symmetry. The main result of our 
calculations is that the properties of the accretion flow do not depend as much on the outer 
boundary conditions (i.e., the amount as well as distribution of the angular momentum) as 
on the geometry of the non- accreting matter. Additionally, we find that the mass accretion 
rate vs. angular momentum relation for a given angular distribution of I has three regimes: 

(i) for low I (i.e., I < IRgc for all 9 at large radii), a torus does not form and Ma — ^b, 

(ii) for intermediate Z, or more appropriately where there is a narrow range of 9 for which 

I > 2Rsc so that 9o > 9t, Ma ~ const and (iii) for high I, or in the case for which I > 2Rsc 
at nearly all 9 so that 9o < 9t, Ma decreases with increasing Iq. Our limited data suggest 
that in the second regime, the actual value of the constant is dependent on the ratio 
Rs/Rb for large Rs/Rb but becomes independent on the ratio for small Rs/Rb (i-e., for 
Rs/Rb ~ 10~^). We conclude that the inclusion of even slow rotational motion of the 
inviscid fiow at large radii can significantly reduce Ma compared to the Bondi rate, as the 
Ma vs. I relation in the third regime indicates. For Rs/Rb > 10~^'^, our results show that 
Ma can be reduced by ~ 1.5 orders of magnitude compared to the Bondi rate. To reduce 
the mass accretion rate more, the accretion funnel needs to be much narrower than the 
funnels allowed by our assumption of low angular momentum (eq. 27) for Rs/Rb > 10^^'^. 
It remains to be seen whether simulations for values of Rs/Rb as low as those observed in 
astrophysical systems (i.e., Rs/Rb ~ 10~^) will confirm our predictions that the accretion 
funnel can be very narrow and subsequently the mass accretion rate very small. 

We note that the discontinuous change from quasi-spherical to a disk-like fiow was 
first found by Abramowicz & Zurek (1981). Abramowicz & Zurek studied analytically 
the adiabatic accretion of a radial flow onto a black hole using the PW potential. They 
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considered the case of constant specific angular momentum and also found that for small I, 
the flow is quasi-spherical and becomes transonic at x^^ ^ R'g while for sufficiently large I, 
the flow has a disk-like pattern and Xg < 3R'g. Applying Abramowicz & Zurek's solution 
for the sonic point, we flnd very good agreement between the location of the sonic point on 
the equator and their prediction. In particular, forn = 1/(7 — 1) = 3/2, I = 2Rsc, and 
zero total energy, eq. 3.5 in Abramowicz & Zurek (1981) yields that the sonic point equals 
1 -|- \/3, which is in a good agreement with our result for 9 — 90° where the sonic radius is 
2.5Rs (e.g., see Figure 4). Our simulations also agree with Abramowicz & Zurek's result 
that the transition between quasi-spherical and disk-like accretion occurs at I ~ 2Rsc. 

We flnd that the shape of the polar funnel through which matter is accreted is partially 
constrained by the torus, which is made of matter that cannot accrete. The torus consists 
of material with roughly uniform angular momentum and its structure does not depend 
much on the outer conditions. The accretion funnel is therefore also not very sensitive to 
the outer conditions in the hmit that it is broad at large radii (i.e., large 60). This is the 
key reason for the insensitivity of Ma to / in the second regime and its weak sensitivity in 
the third regime: the mass accretion rate depends on the geometry of the sonic surface 
at radii where the presence of the torus is important, and not only on the geometry at 
the outer boundary. Thus, the accretion flow in the funnel should depend on the physics 
that controls the flow in the torus. In particular, the introduction of energy dissipation, 
and the transport of energy and angular momentum in the torus, may change the shape 
of the torus and its effect on the polar regions, where material can accrete without the 
transport of angular momentum. One could argue that the total mass accretion rate (via 
both the funnel and torus) onto the black hole should increase when transport of angular 
momentum and consequently accretion via the torus are allowed. However, it is not clear 
by how much Ma will increase, if at all, because the energy and angular momentum from 
the torus are likely to be deposited in the polar region (sec e.g., Stone, Pringlc & Bcgclman 
1999; Blandford & Begelman 1999; Blandford & Begelman 2002a, 2002b; Hawley & Balbus 
2002). Indeed, outfiow from the torus may interfere with the infiow in the funnel and M^ 
may well decrease. There is also a possibility that accretion via the funnel will decrease, not 
because of the dynamical effects due to outfiowing material but because a torus in which 
energy is dissipated should be hotter and thicker than a torus in which energy cannot be 
dissipated. On the other hand, it is not clear how accretion via the torus will change when 
supersonic accretion in the funnel is present. Our simulations show that the range of 6 that 
is occupied by the accretion funnel increases with decreasing distance from the black hole. 
In particular, accretion through the inner radius occurs for the entire range of 6. We expect 
that the presence of a supersonic accretion fiow near the equator at small radii may cause 
accretion via the torus to be less effective compared to the situation in which the material 
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in the polar region is static. Most simulations of the formation of the torus have assumed 
that the material near the poles, outside the initially hydrostatic rotating torus, is static or 
subsonic, and therefore unimportant dynamically. 

Even within our simple framework of inviscid flow, we anticipate that our results may 
change if we relax some of our assumptions about the geometry at the outer boundary. 
The assumptions we adopted are as simple as possible because we do not know the 
magnitude of the angular momentum of material in the environment of SMBHs, nor the 
angular distribution of /. In particular, it is rather unlikely that the / distribution is 
axisymmetric. Therefore, fully 3D calculations are required to explore how Ma responds 
to non-axisymmetric / distributions. In the context of inviscid HD calculations similar 
to ours, one would expect that 3D effects may reduce Ma compared to 2D axisymmetric 
calculations. For example, it is possible that rotation at large radii occurs not just around 
one axis but around two or more axes. It is also possible that rotation occurs around an 
axis that changes with r. Additionally it is plausible that there are "winds" flowing past the 
SMBH or that the entropy in the SMBH environment is variable. In such a case, material 
with too high / to be accreted may occupy nearly the entire range of 9 at small radii (near 
or inside the sonic radius) and prevent low-/ material from accreting. In terms of accretion 
funnels, this would correspond to the situation in which the 'funnel' is very narrow and 
maybe not lined up with a radial vector. 

Perhaps the best studied massive black hole with a very low luminosity is the black 
hole in the Galactic center. Many models and ideas for how to explain very low SMBH 
luminosities have been explored in the context of Sgr A*. In particular, Melia (1992; 
1994) proposed a spherical accretion model in which the accretion flow is assumed to be 
in free-fall until a Keplerian disk is formed within a small circularization radius. Coker & 
Melia (2000) looked at the problem of spherical accretion but with the magnetic field being 
in subequipartition. The latter results in a reduced bremsstrahlung emissivity and can help 
to explain the low luminosity of Sgr A*. A relatively small distance to Sgr A* allows us 
to map the vicinity of the black hole at the Galactic center. The complexity of such maps 
(e.g., in radio) motivate three dimensional HD simulations. For example, Coker & Melia 
(1997) performed three-dimensional simulations of Bondi-Hoyle accretion of stellar winds 
onto a black hole. Glearly Sgr A* shows us that accretion onto black holes is a complex 
phenomenon and the Bondi accretion formula should be used with great caution because 
the assumption of spherical accretion is most certainly an oversimplification. 

We finish with the observation that after three decades of studying the Bondi problem 
with shght rotation at large radii (see, e.g., Henriksen &: Heaton 1975, Lynden-Bell 1978; 
Cassen & Pettibone 1976; Sparke & Shu 1980; Sparke 1982; Abramowicz & Zurek 1981, 
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for analytic attempts to solve this problem and see references in Section 1 for 2D and 3D 
numerical simulations), we still find new complexities in the behavior of accretion flows. 
In our next phase of studying these flows, we will consider 3D MHD models of radiatively 
inefficient flows. 
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Fig. 1. A sequence of specific angular momentum contours (top) and velocity fields 
(bottom) from run B04fla at times 0.03, 0.06, 4.56, and 9.27. The specific angular 
momentum is in units of 2RsC. The minimum of / (contour closest to the rotational 2;— axis) 
is 0.3 and the contours levels are equally spaced at intervals of Al = 0.3. Note that the 
maximum of I is 3.6 at the beginning of the simulation and 3.3 at the end of the simulation, 
indicating that lower angular momentum material is displacing higher angular momentum 
material in the torus. Only the I — 0.9 contour is labeled to show the boundary between 
the flow that does and does not accrete. We suppressed the velocity in the innermost part 
of the flow to better show the flow pattern at large radii. 

Fig. 2. The specific angular momentum contours and velocity field at the end of 
run B04fla. This figure is an enlargement of the rightmost panels from Figure 1. Note 
that at small radii the nonaccreting flow, the torus, is of nearly constant specific angular 
momentum. The "figure- 8" contour on both panels marks the radial sonic surface (the 
location where Mr — Vr/cg — 1). 

Fig. 3. Quantities at the inner boundary in model B04fla at tf — 9.27. The 
solid and dashed lines on the top panel show the radial and latitudinal Mach numbers, 
respectively. The specific angular momentum (middle panel) is in units of 2Rsc while 
the mass flux density (bottom panel) is in units of the Bondi mass flux density [i.e., 
-{pVr)B = MB/{Anr^)]. 

Fig. 4. Quantities on the sonic surface in model B04fla at tf — 9.27. The top panel 
shows the sonic radius as a function of the latitude in units oi Rb- The second panel from 
the top shows / in units of 2Rsc (solid line) and the latitudinal Mach number (dashed line). 
The second panel from the bottom shows the mass flux density in units of the Bondi mass 
flux density at the radial sonic point, [i.e., — (pfr)_B = ^B/(47ri?^a;^)]. The bottom panel 
shows a measure of the deviation from radial flow, A^, as deflned in the text. 

Fig. 5. The time evolution of the mass accretion rate in units of the Bondi rate, for 
model B04fla. 

Fig. 6. The specific angular momentum contours at the end of run BOlfla. The 
right-hand panel is an enlargement of the left-hand panel. The "figure-8" contour on the 
right-hand panel marks the the radial sonic surface (the location where Mr = 1). The 
spacing of the / contours is as in Figure 1. Note that runs B04fla (Fig. 1) and BOlfla 
have different / distributions at large radii yet at small radii the ffows in both runs are 
qualitatively similar (compare figure 2 and the right-hand panel here). 

Fig. 7. The mass accretion rate as a function of 60- The angle 60 is the polar angle 
at which / = 2Rsc at the outer boundary. The solid line represents results for all of our 
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Table 1: Summary of parameter survey. 
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B Comments 



To = 12Rb 



* We use the following convention to label our runs: the first character in the name refers to R'g, (i.e., A, 
B, and C are for R'g = 10-2, 10^^^ and 10"3-5^ respectively). The second and third characters refer to 6o 
(i.e., 01, 02, 03, 04, 05, 06, 07, 08, 09 and, 10 are for Oo = 90°, 68°, 53°, 44°, 29°, 24°, 18°, 14°, 8.1°, 6.3°, and 
4.5°, respectively), the fourth and fifth characters refer to the angular distribution of I (i.e., /I, /2, and /3 
are for 1 — | cos 6*1, 1 — cosi° 6*, and the step function (see eq. 25), respectively), and finally the sixth 



character refers to the outer radius of the computational domain (i.e., a and b stands for Tq 

12 Rb, respectively). 



1.2 Rb and 
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model with Rs/Rb — 10~^, regardless of the angular distribution of I and regardless of 
Iq (see Table 1). The dashed line represent the theoretical prediction that scale with 
the solid angle within which I < 2Rsc at Rb- The mass accretion rate is in units of the 
corresponding Bondi rate. 
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